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Abstract 

We apply the finite size scaling analysis to the derivative of the density of the effec- 
tive action for the lattice U(l) pure gauge theory in an external constant magnetic 
field. We found the presence of a continuous phase transition. Moreover, our exti- 
mate of of the critical parameters gives values consistent with those extracted from 
the analysis of the specific heat. 



1 Introduction 



In a recent paper [1] (henceforth referred as I) we proposed a new method 
that allows to evaluate the gauge invariant effective action for external back- 
ground fields in lattice gauge theories. In I we considered the so-called lattice 
Schrodinger functional [2]: 



Z 



A 



cxt 



(A ext \exp(-HT)V\A cxt ) = JvUexp(-S w ) 



where Sw is the standard Wilson action. In Equation (1) we integrate over 



the links U^{x) with the constraints 



u,(x)\ X4=0 = u; x \ (2) 

where 

U e *\x) = P exp l^+iag J dt A)f (x + atjl)\ . (3) 
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The lattice effective action for the background field A cxt (x) is given by: 



cxtl 



icxt 



In 



Z\U_ 



(4) 



where T is the extension in the Euclidean time. In the continuum where 



T — > oo it is easy to show that T A 
presence of the background field A ext 



ext. 



reduces to the vacuum energy in 



x). Note that our effective action is by 
definition gauge invariant. Moreover the effective action (4) can be used for a 
non-perturbative investigation of the properties of the quantum vacuum. 



2 U(l) in a constant background field 



In I we considered the four-dimensional U(l) lattice gauge theory in an exter- 
nal constant magnetic field: 



AJ?ffl = 6 k , 2 x 1 B. 

The lattice links corresponding to the vector potential (5) are 

£7| xt (a;) = exp [iagBxi] = cos(agBxi) + i sin(agBxi) , 
U^\x) = U™\x) = UZ\x) = 1 , 

where, due to the periodic boundary conditions, we have 
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a z gB 
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-n 



cxt 



(5) 



(6) 



(7) 



with n ext an integer, and L\ the lattice extension in the x\ direction (in lattice 
units). 

Due to the gauge invariance the effective action only depends on the constant 
field strength Fffi. Thus we consider the density of the effective action 



A 



ext. 



A 



cxt 



V 



In 



cxtl 



Z(0) 



(8) 



where Q — V ■ T, with V the lattice spatial volume. In order to avoid the 
problem of computing a partition function, we focused on the derivative of 



text 



with respect to (3: 



text 



- cos9^(x)\ -(^2 cos^(x)\ 



(9) 
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From our previous discussion it is clear that e 



icxt 



is the vacuum energy 



density. Thus we expect that e' 



A c 



behaves like a specific heat. Indeed, in 



A 



ext, 



displays a maximum 



I we found that, near the critical region (3 pa 1, e' 
which increases by increasing the lattice volume. In addition, the peak shrinks 
and shift towards (5 — 1 by increasing L = Q 1 / 4 . 



3 Finite Size Scaling Analysis 



The aim of this paper is to determine the critical parameters by applying the 
standard finite size scaling analysis [3]. Consider a finite system with effective 
linear dimension L. In general the thermodynamical quantities depend on 
£/a and L/a, where £ is the correlation length and a is the lattice spacing. 
The finite size scaling hypothesis assumes that near the critical point the 
microscopic length a drops out. 

Let us consider, for instance, the generalized susceptibility x- I* 1 a infinite 
system near the critical point we have 

x~r /v - (io) 

In the finite geometry with size L, according to the finite size scaling hypoth- 
esis, we can write 

X = C /u ^/L) (11) 
where the scaling function (p{u) is smooth. Now, from £ ~ \T — T c \~ v it follows 

X = L^^(L^(T-T C )) , (12) 

where (f>(u) is another smooth scaling function. On a finite system x is sup- 
posed to be analytic in T. Thus <f>(u) must be a smooth function with a peak 
at some value u = uq. Near the peak we parametrize <j>(u) as: 

= 7 % 7 i ( 13 ) 
a 2 {u - u Q y + 1 

with ai, a 2 , u positive constants. 

From Equations (12) and (13) we obtain easily: 

= aiL^ 
X a 2 L 2 /»[T -T*(L)] 2 + 1 1 J 
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where 



T*(L) 



-1/1/ 



(15) 



Then, if we plot the susceptibility versus the temperature, we see a peak at 
T*(L) with a width L^ 1 ^ and height ~ L 1 ^ . Further, the critical exponents 
7 and nu are related by the so-called hyperscaling [3] relation [3] 



7 = 2 - vd 

where d is the dimensionality of the system. Equation (16) implies that 



1 

v 



d 
2 



2v 



(16) 



(17) 



These considerations can be directly applied to our generalized susceptibil- 



ity e' 



A 



ext 



In order to determine the critical parameters we need to simu- 



late our gauge system on lattices with different sizes. However we must take 
fixed the external magnetic field, for e' A cxt depends even on gB. Thus, 
according to Eq. (7), we must work with L\ and n cxt fixed. As a conse- 
quence we perform the numerical simulations on lattices with sizes L\ = 64, 
L 2 = L 3 = L 4 = L = 6, 8, 10, 12 J 14 and n cxt = 2. It is worthwhile to stress 
that the contributions to e 



l j^ext 



due to the frozen time slice at x 4 = and 
to the fixed boundary conditions at the lattice spatial boundaries must be 
subtracted. Indeed, the physically relevant quantity is due to the "internal" 
links, i.e. only the dynamical links must be taken into account in evaluating 



A c 



In Figure 1 we display the derivative of the internal energy e' 



A c 



versus (3 for the 64 x 12 3 lattice. Note that in the strong coupling region (3 < 0.7 
A oxt vanishes so that the density of the effective action can be recovered 
by integrating over (3 without any subtraction. 



In Figure 2 we show e' A ext 
sizes. We try a fit of the pea 

e' \A cxt 



near the critical region for the different lattice 
<; according to Eq. (13): 



ai(Asff) 



a 2 (L cS )i(3-(3*(L cS )]z + l 
where the effective size of the lattice is [3] 



n 



1/4 
int 



(18) 



(19) 



In Equation (19) fli nt is the internal lattice volume, i.e. the volume occupied 
by the dynamical lattice sites. If we denote by f2 cxt the lattice sites whose links 
are fixed according to Eq. (6), thus we get: 



Q cxt = L X L 2 L Z + (L 4 - \){L X L 2 L Z - {L x - 2){L 2 - 2)(L 3 - 2)) , 



(20) 
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where the first term is the volume of the frozen time slice X4 = while the sec- 
ond term is due to the fixed links at the lattice spatial boundaries. Whereupon 
we have: 



int 



-^1-^2-^3-^4 — ^cxt — ^ — ^ext 



(21) 



It turns out that our lattice data are in satisfying agreement with Eq. (18). 
Indeed near the pseudocritical coupling we fit the data according to Eq. (18) 
by means of the standard MINUIT subroutine of the CERN library. We find 
rather good fits (solid lines in Fig. 2) with x 2 /d.o.f. ~ 1. As a result of the fits 
we get the values of the three parameters for the five different values of L c g. 
According to the finite size scaling hypothesis ai(L e fr) should behave as 



ai(L efT ) = cLl£ . 



(22) 



In Figure 3 we report ai(L e ff) versus L e g together with the result of the fit 
Eq. (22). We find 



c= 0.63(24) 
= 0.84(14) 



(23) 
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with x 2 /d.o.f. 
Eq. (15): 



1.2. The pseudocritical couplings are fitted according to 
P*(L cS )=[3 c + dL- s 1/u . (24) 



We find (see Fig. 4): 



v 



3.48(65) 



(25) 



(5 C = 1.0089(47) 



with % 2 /d.o.f. = 1.7. The values Eqs. (23) and (25) for the critical exponent are 
quite consistent with the hyperscaling relation Eq. (16). Indeed from Eqs. (17) 
and (23) we find 



- = 2.42(7) , 



(26) 



which is compatible with Eq. (25) within two standard deviations. 
Finally, in Fig. 5 we check that the universality law in the critical region: 



-^eff £ int 



icxt 



I%{P-P C ) 



(27) 



holds quite well for the lattices with L = 10, 12, 14 corresponding to L eff = 
13.75, 16.16, 18.46 respectively. 
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4 Conclusions 

We would like to emphasize that our extimation of the infinite volume critical 
coupling is in perfect agreement with the values extracted from the specific 
heat [4,5]. However it is well established that on lattices with periodic bound- 
ary conditions, the U(l) phase transition is exhibiting a weakly first-order 
nature [5]. On the other hand, our data are consistent with a continuous 
phase transition. As a matter of fact, it turns out that the phase transition 
of compact U(l) pure lattice gauge theory with the Wilson action on lattices 
with closed topology [6,7] or with fixed boundary conditions [9] is continuous. 
Moreover our extimation of the critical exponents is consistent with the results 
of Ref . [6] . Thus, the finite size scaling analysis discussed in this paper results 
in a successfully quantitative test of our proposal of a gauge invariant effective 
action on the lattice. 
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FIGURE CAPTIONS 



Figure 1. The derivative of the energy density due to links belonging to f2 int versus 

P for the 64 x 12 3 lattice. 
Figure 2. e' int near the critical region for the lattices 64 x 6 3 (open circles), 64 x 8 3 

(open squares), 64 x 10 3 (open triangles), 64 x 12 3 (open diamonds), and 

64 x 14 3 (open crosses), with the superimposed fits Eq. (18). 
Figure 3. The maximum of e[ nt obtained from the fit Eq. (18) versus L c r with 

superimposed the fit Eq. (22). 
Figure 4. The pseudocritical couplings at various lattice sizes with superimposed 

the fit Eq. (24). 

Figure 5. The universality law Eq. (27) in the critical region. Symbols as in Fig. 2. 
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Fig. 1. The derivative of the energy density due to links belonging to Q- m t versus (3 
for the 64 x 12 3 lattice. 
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Fig. 2. e[ nt near the critical region for the lattices 64 x 6 3 (open circles), 64 : 
(open squares), 64 x 10 3 (open triangles), 64 x 12 3 (open diamonds), and 64 x 
(open crosses), with the superimposed fits Eq. (18). 
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P 



(3 (L eff )= p c + d L 



•1/v 
eff 



1.25 



1/v = 3.48(65) 
p c = 1.0089(47) 
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Fig. 4. The pseudo-critical couplings at various lattice sizes with superimposed the 
fit Eq. (24). 
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5. The universality law Eq. (27) in the critical region. Symbols as in Fi 
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